Subband filling and Mott transition in Ca 2 _ x Sr x Ru04 
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A new concept is proposed for the paramagnetic metal insulator transition in the layer perovskite 
Ca2- a; Sr 2 :Ru04. Whereas the pure Sr compound is metallic up to very large Coulomb energies 
due to strong orbital fluctuations, structural changes induced by doping with Ca give rise to a 
interorbital charge transfer which makes the material extremely sensitive to local correlations. 
Using dynamical mean field theory based on finite temperature multi-band exact diagonalization 
it is shown that the combination of crystal field splitting and onsite Coulomb interactions leads to 
complete filling of the d xy band and to a Mott transition in the half-filled d xz , yz bands. 
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The layer perovskite C^-aSr^RuC^ has attracted 
wide interest during recent years because of the complex 
sequence of electronic and magnetic phases which arise 
when Sr is iso-electronically substituted by Ca 0, 0, Q . 
While the pure Sr compound exhibits unconventional su- 
perconductivity Q , the structural distortions induced by 
the smaller Ca ions ultimately lead to an antiferromag- 
netic Mott insulator. The physical mechanism of how 
this multi-band material with four electrons per Ru ion 
evolves from the metal Sr2Ru04 towards the insulator 
Ca2Ru04 is presently not well understood. 

According to band structure calculations [f| the t 2g 
orbitals in Sr2Ru04 are approximately equally occupied. 
Because of the planar geometry, these bands split into a 
wide, nearly two-dimensional d xy band and two narrow, 
nearly one-dimensional d XZ}VZ bands. Clearly, the Mott 
transition in such a highly anisotropic system cannot be 
understood in terms of single-band models, or multi-band 
extensions assuming degenerate subbands. To analyze 
the role of local Coulomb interactions in Ca2- a; Sr :E Ru04 
several aspects were investigated within dynamical mean 
field theory (DMFT) Q, using a variety of impurity 
solvers to treat onsite correlations 0, @, & 0, EH E3] • 

While it is generally agreed upon that Coulomb inter- 
actions in Sr2Ru04 lead to a sizeable band narrowing and 
effective mass enhancement, and a shift of the d xy van 
Hove singularity towards the Fermi level due to d XZ}V 



band model consisting of wide and narrow semi-circular 
densities of states. To account for the Ca induced octa- 
hedral distortions the d xy band was assumed to be 
narrower than the d XZtVZ bands. Accordingly, in contrast 



to Ref. 



d xy is the first band to become insulating 



with increasing U. An important parameter in this work 
is the crystal field splitting A between t 2g orbitals, which 
yields Mott transitions at much smaller values of U than 
for A = 0. In both models, a second interorbital charge 
transfer due to additional structural modifications, i.e., 
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(0.5,0.75,0.75) -> (1.0,0.5,0.5), 



quired in the limit x — > to yield the antiferromagnetic 
insulating properties of Ca2Ru04 [8j|. 

In the present work, we use finite temperature exact 
diagonalization (ED) DMFT [l5[ to study the nature of 
the Mott transition in Ca2- a; Sr : i;Ru04. Since a t 2g tight- 
binding Hamiltonian including the full complexity of the 
octahedral distortions as a function of Ca/Sr concentra- 
tion is not yet available, we use the Sr2Ru04 density of 
states components as a single-particle starting point (see 
Fig. 1). As shown by Fang et al. 14|, a key effect caused 
by rotation, tilting and flattening of oxygen octahedra is 
the increasing d xy orbital occupancy with increasing Ca 
concentration: In the limit x = (pure Ca), the spin- 
averaged occupancy is about n xy sa 0.83, compared to 
n xy ks 0.64 for x — 2 (pure Sr). To account for this 
charge transfer, we allow, as in Ref. [l2j], for a crystal 



d xy charge transfer @, d 0, EH , the Mott transition near field splitting A between d xy and d xz>yz states. For in- 



*xy 

the Ca end of the phase diagram is complicated because 
of the charge rearrangement among t 2g orbitals when the 
Ca/Sr concentration is varied. Using the non-crossing 
approximation fl3| . Anisimov et al. [<§] obtained succes- 
sive, 'orbital selective' Mott transitions upon increasing 
the onsite Coulomb energy U: first for the narrow d xz , yz 
bands and subsequently for the wide d xy band. These 



stance, a lowering of the d xy bands by A = 0.2 (0.4) eV 
yields n xy « 0.74 (0.83). The total occupancy of the t 2g 
bands is four, independently of Ca/Sr concentration. 

The main result of this work is a new mechanism for 
the Mott transition in this multi-band system: For re- 
alistic values of the Coulomb energy U and crystal field 
splitting A, we find a gradual filling of the d xy subband 



transitions arise via a d 



ij^xy i 



xy 



*>x Z ,yz 



charge transfer, i.e., with increasing U. Once the d. 
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) « (2/3,2/3,2/3) -> (0.5,0.75,0.75). Or- is complete, i.e., (n 
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xz,y Z d xy charge transfer 
n yz ) (1.0,0.5,0.5), a Mott 



bital selective Mott transitions, with the same interor- 
bital charge transfer, were found also by Dai et al. [T3| 
within slave boson mean field calculations for a three- 



transition takes place in the remaining half-filled d xz _ yz 
bands. The critical Coulomb energy at which this tran- 
sition occurs is lower for larger crystal field A. 
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FIG. 1: Upper panel: 
Solid (red) curve: wide 



Density of states of S^RuCU. 
d x 



band; dashed (blue) curve: 
narrow d xz , yz bands (Ef = 0). Lower panel: Fermi 
surface (schematic) of Sr2Ru04 (dashed blue curve) and 
Cai.5Sro.5Ru04 (solid red curve). In the former (latter) case, 
the d X y van Hove singularity at M is above (below) Ef- 



This scenario is differs qualitatively from the ones pro- 
posed in Refs. (slH^|: Instead of successive transitions in 
the half-filled d xy band and in the 3/4-fflled d XZiVZ bands, 
we find a common d xy band filling and d XZiVZ Mott tran- 
sition. The d xy band filling, induced by the combined 
effect of crystal field and Coulomb interactions, signifies 
a change from a correlated metallic state to a correlated 
band insulating state and does not correspond to a Mott 
transition. On the other hand, since a full d xy band im- 
plies half-filled d xz>yz bands, the latter readily undergo 
a standard metal insulator transition. Thus, the present 
three-band system does not exhibit orbital selective Mott 
transitions. Our solution solely relies on the structure 
induced increase in d xy occupancy 14| and is consistent 



with the (n xy , n. 



) = (1.0, 0.5, 0.5) configuration ob- 



| . Thus, we do not 
2/3 -> 1/2 -> 1 as 



tained for Ca2Ru04 within LDA+U 
require the more complex path n xy r 
Sr is replaced by Ca. 

The details of our multi-band ED approach are given 
in Ref. [l5j |. Exploiting the sparseness of the impurity 
Hamiltonian, much larger clusters representing the bath 
can now be handled, thereby significantly reducing finite 
size effects. All calculations are carried out for full Hund 
exchange J. To obtain the orbital occupancy n, as a 




FIG. 2: Orbital occupancy rij as a function of Coulomb 
energy for various crystal field splittings A, derived within 
ED/DMFT for T = 20 meV. Solid (red) dots: n xy , empty 
(blue) dots: n XZtyz . The lines are guides to the eye. 



function of U, for a particular A, we assume J = [7/4. 
The interorbital Coulomb energy is U' = U — 2J. To 
achieve fast convergence we use the cluster size n s = 9 
(2 bath levels per impurity orbital); for greater precision 
specific points are calculated with n s = 12. The temper- 
ature is assumed to be T = 20 meV. We focus here on 
the paramagnetic Mott transition. 

Fig. 1 shows the ti g density of states components for 
Sr2Ru04, derived from a tight-binding Hamiltonian 0] 
fitted to the LDA band structure of Ref. [j| . The band 
extrema were adjusted slightly in order to accomodate 
four electrons below Ep. The d xy van Hove singularity 
lies about 60 meV above Ep. The lower panel shows the 
Fermi surface. The 7 sheet corresponds to d xy , the a and 
j3 sheets to d xz>yz . The latter do not cross because of a 
small hybridization term. Also shown is the Fermi surface 
for a crystal field A = 0.1 eV. Since the d xy van Hove 
singularity now lies below Ep, the 7 sheet has turned 
from electron-like to hole-like. 

Fig. 2 summarizes the variation of the ti g orbital oc- 
cupancies with onsite Coulomb energy U. In the case of 
Sr 2 Ru04 (A = 0), a gradual charge transfer from d xz ,yz 
to d xy states is found, in agreement with the trend found 
within QMC/DMFT for small U [7]. The system is seen 
to remain metallic up to very large U. This result is con- 
sistent with previous calculations 0] which revealed no 
Mott transition up to U = 4 eV. Near U c ~ 10 eV the in- 
terorbital charge transfer is complete with n xy = 1. The 
remaining half-filled d xz ^ yz bands then undergo a Mott 
transition (see below). 

Note that the filling of the d xy band for A = accel- 
erates towards increasing U. The reason for this trend 
is that the larger d xy occupancy reduces the intra-£2g 
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screening of U. Since at the same time the occupancy 
of the d xz ,yz bands approaches one-half, the tendency for 
the latter bands to undergo a Mott transition becomes 
rapidly more favorable. In the absence of the d xy band, 
the half-filled d XZtVZ bands (width W — 1.2 eV) would ex- 
hibit a common metal insulator transition at U C — 1A eV, 
i.e., far below U = 10 cV (U c w 0.8 W for twofold de- 
generate semicircular bands and J = U/4 Thus 
the results for A = may be interpreted as d xz ^ yz Mott 
transition delayed by strong orbital fluctuations. 

In view of this picture it is plausible that, at finite Ca 
concentrations, the greater initial occupancy of the d xy 
band associated with A > gives rise to complete d xy 
filling and a d xz ^ yz Mott transition at progressively lower 
values of U c . This is confirmed by the results shown in 
Fig. 2. According to Fang et al. [14[ the crystal field split- 
ting for Ca2Ru04 induced by octahedral distortions is 
about 0.4 eV. The d xy band filling and d XZtVZ Mott tran- 
sition is then shifted to U c ~ 4 eV. Thus, the reduced or- 
bital fluctuations greatly diminish the delay of the Mott 
transition in the half-filled d xz>yz bands. A slightly lower 
U c would be obtained if other effects neglected so far, 
such as distortion induced band narrowing and interor- 
bital hybridization [l3 |. are taken into account. 

In the range 0.5 < x < 2.0 the octahedral distortions 
consist mainly of rotations about the z-axis. The as- 
sociated lowering of the d xy band is rather small, giv- 



ing A < 0.2 eV |14|. According to the results shown 
in Fig. 2, the correlation induced d xy filling and d xz%yz 
Mott transition then occur at U c « 6 eV. In this doping 
region the system therefore remains metallic. This result 
is consistent with angle resolved photoemission data 
and optical data 0] which show that near x = 0.5 all ti g 
states are itinerant. 

A crystal field A = 0.1 . . . 0.2 eV is large enough to 
push the d xy van Hove singularity below Ep. This also 
agrees with the photoemission results Q which yield 
nearly the same Fermi surface for Cai.5Sr .5RuO4 and 
Sr 2 Ru04, except that, because of the lowering of the van 
Hove singularity, the 7 sheet has changed from electron- 
like to hole-like, as indicated in Fig. 1. These data 
are naturally explained by the results shown in Fig. 2, 
whereas they are difficult to reconcile with the trend to- 
wards half-filled d xy bands proposed in Refs. (j. 12|. 



To prove that d xy band filling and d xz ^ yz Mott tran- 
sition coincide we show in Fig. 3 the ti g quasi-particle 
spectra for A = 0.4 eV. To avoid uncertainties stem- 
ming from the extrapolation from Matsubara frequencies 
to real frequencies we give here the spectra of the clus- 
ter Green's functions, Ai(u>) = — -ImGi(w + iS), with 
6 = 50 meV. While spectral details differ from those of 
the solid [l5| > the cluster results are adequate for the dis- 
tinction between metallic and insulating behavior. The 
spectra for U = 3.0 eV reveal that all subbands are 
metallic and exhibit appreciable spectral weight below 
the single particle bands, associated with lower Hub- 
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FIG. 3: Quasi-particle spectra calculated using ED/DMFT 
for crystal field splitting A = 0.4 eV; T = 20 meV. (a) metal- 
lic region U = 3.0 eV; (b) transition near U c ps 4.2 eV; (c) 
insulating region U = 4.5 eV. Solid (red) curves: d xy band, 
dashed (blue) curves: d XZtVZ bands. 



bard bands. According to Fig. 2, the occupancies are 
n xy — 0.84 and n xz<yz = 0.58. In contrast, at U = 4.5 eV 
the d xy band is filled and the half-filled d XZiyz bands ex- 
hibit a clear separation into upper and lower Hubbard 
bands. The transition between these two regions occurs 
at U c w 4.2 eV. Note that the insulating gap arises be- 
tween the filled d xy band and the upper Hubbard bands 
of the half-filled d xz%yz bands. The quasi-particle proper- 
ties near this transition will be discussed elsewhere. 

We point out that, in general, the d xy band filling and 



M xz,yz 



Mott transition do not need to take place at the 
same U. For instance, if the d XZjVZ bands were much 
wider, their Mott transition would occur above the d xy 
band filling. Also, a stronger crystal field could give a d xy 
filling at very small U, and a Mott transition in the half- 
filled d xz . yz bands at a larger U. Because of the small 
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width of these bands the trend seen in Fig. 2 suggests 
that this possibility should arise only for A > 0.5 eV. An 
analogous effect was discussed by Manini et al. [T7| for a 
two-band model with equal bands of semi-circular density 
of states, off-set via a crystal field A. For unit total 
occupancy and small A, a single transition was found 
where one band is pushed above Ep and the other (then 
half-filled) exhibits a Mott transition. For larger A, one 
band is emptied at small U while the second remains 
metallic up to a Mott transition at larger U. 

We also note that the Mott transitions in the 3d 1 
perovskites LaTiC>3 and YTiC>3 studied by Pavarini et 
al. [H} reveal an almost complete emptying of two t2 g 
subbands, and a metal insulator transition in the re- 
maining nearly half-filled band. Evidently, in all of these 
multi-band systems, the Mott transition is made feasible 
by a striking suppression of orbital fluctuations. 

In the work discussed above the d X z,yz bands are 
degenerate. In real Sr2Ru04 these orbitals interact 
weakly, giving slightly different subband densities of 
states N xz ± yz (uj) of identical width. Ca induced dis- 
tortions will enhance these differences, so that the Mott 
transition in these subbands becomes non-trivial. This 
will be also addressed in future studies. 

To analyze the Mott transition in Ca2- :c Sr :c Ru04 we 
have focused on the variation of the subband occupancies 
with Coulomb energy. In reality, U should be roughly 
constant as a function of x, with U w 3.1 eV and 
J w 0.7 eV according to constrained LDA calculations for 
x = 2 [11]. Thus, in Fig. 2 a vertical line near U — 3 eV 
should qualitatively cover the low temperature phase dia- 
gram. A = with n xy — 0.70, n XZy . 



0.65 corresponds 



to metallic Sr2Ru04. A = 0.2 eV with n x 



0.76, 



terial becomes highly sensitive to local correlations. In 
the Ca rich compound orbital fluctuations are sufficiently 
weak that the Mott transition occurs at realistic values 
of U. On the basis of this picture it would be very inter- 
esting to study the region x < 0.5 more closely in order 
to understand the orbital selective mass enhancement Q 
and the rich magnetic phases of this material [l[. 
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n xz ,yz — 0.62 represents the case Cai.sSro.5Ru04, which 
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